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ABSTRACT 

It is shown that the renormalisation group flow in coupling constant space can be 
interpreted in terms of a dynamical equation for the couplings analogous to viscous 
fluid flow under the action of a potential. For free scalar field theory the flow is 
geodesic in two dimensions, while for D 7^ 2 it is only geodesic in certain limits, e.g. 
for vanishing external source. For the 1-D Ising model the renormalisation flow is 
geodesic when the external magnetic field vanishes. 
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Introduction 



The machinery of the renormahsation group has evolved steadily over the 40 years 
since its inception and recently the notion of a geometry on the space of couplings has 
received some attention [1] [2] , (for a review in a thermodynamic context see [3] ). 
The aim of this paper is to use these geometrical notions to interpret the /3-functions 
as a dynamical flow and search for a dynamical equation governing the evolution of 
the couplings of a quantum field theory under changes of scale. It will be argued that, 
with a specific choice of metric on the space of couplings, the fiow of the couplings is 
analogous to fluid flow in a curved space, under the influence of a potential, including 
friction. Under certain circumstances the flow is geodesic, with the constraint that 
the kinetic energy equals the potential energy. 

For example for free massive fleld theories coupled to a constant source, the mass 
and the source j parameterise a two dimensional space and it will be shown that 
the renormahsation group flow is geodesic along the j-axis (m^ = 0) and along the 
m^-axis (j = 0), but not otherwise, (except in two dimensions which is special in that 
all trajectories are geodesies). 

The concept of geodesic flow requires introducing a metric on the space of couplings 

[2] [4] [5] . The significance of a metric was highlighted by Zamolodchikov who utilised 
the concept in the proof of the c-theorem in two dimensional field theory [4] . A related 
metric in D-dimensions was considered by O'Connor and Stephens in [2] and this is 
the metric which will be studied here. 



The construction used by O'Connor and Stephens was to consider the partition func- 
tion Z and free energy W to be differentiable functions on the n-dimensional space 
of couplings, parameterised by g"' ,a — 1, . . . , n. 
Then from the identity 

'"i?V'e-^[^l+^ (1) 



follows the formula 

dW = (dS) (2) 

where dS = (J ^ a{x)d^ x)dg°' is an operator value one - form. O'Connor and Stephen's 
metric is 

{{dS-dW)®{dS-dW)). (3) 

In order to be able to pass to the infinite volume limit densities will be used here. 
Defining 

^a{x) = ^aix) - {^a{x)) , (4) 



we shall take 



to be the metric. 



Gab = J {^a{x)MO))d''x (5) 
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In special co-ordinates in which the action is hnear in the couphngs, Gab can easily 
be derived from the free energy density w, defined via W = J wd^x, from 

Gab = -dadbw. (6) 

Note however that equation (6) is not covariant under general co-ordinate transforma- 
tions, whereas equation (5) is. In general there may be singularities in {^a{x)^b{y)) as 
|a; — 1/| — > which are strong enough to render the integral in (5) divergent and it will 
be assumed that these can be regularised. It will also be assumed that ($a(a;)$6(y)) 
falls of fast enough for large — j/| for (5) to be finite. The formalism will now be 
illustrated with two examples - free field theory and then the 1-D Ising model. 



Free Field Theory 



Consider free field theory in D Euclidean dimensions, with action 

S{cp;j,m^) = yd^a;|^(^(-n2 + m2)(^ + jv|. (7) 



There are two parameters in this theory, j (which shall be taken to be independent 
of position for simplicity) and m?. The partition function is 

Z{j,m^) = / Dipe-^^^' (8) 



Performing the Gaussian functional integral gives the free energy density, using di- 
mensional continuation with dimension D, 



w{j,m ) = - — 



2\ D/2 -p/D D\ 



An J L>(f-1) 2m2" 



(9) 



Since the action is linear in j and the metric as given by (6) is easily calculated 
giving 



ds'^ = -idadbw)dx''dx^ = At df + 



r(2 



2(47r) 



(dm^)^ ^ djdrm? . 



(10) 



As pointed out in [2], it is convenient to use = = —j/m?, rather than j itself 
as a co-ordinate, because this diagonalises the metric. This leads to 



1 



(47r) 



D/2 



D 



{m^)^-' r 2 - — \ {dm^) 



,2\2 



(11) 



To expose the geometry more clearly a further change of co-ordinates is useful. 
Let 




then 



ds^ = dr^ + r^de^. (13) 



This co-ordinate transformation is singular for D — 4, but perfectly regular for 
< D < 4. It is clear from (13) that D = 2 gives the flat metric in r — ^ space, 
but all other dimensions lead to a curved metrics. The curvature scalar is 

which is equivalent to the large volume limit of [2]. 

Note that 7^ = in D = 4 for finite m^, though it is positive for finite r, due to the 
singular nature of the co-ordinate transformation (12) in four dimensions. For D = 0, 
TZ — —(1/2) for finite m^, giving the Lobachevski plane for the Gaussian distribution 
in ordinary statistics, [6] . It is, of course, which is the physical parameter. 

We shall now turn to a discussion of the renormalisation flow in this geometry. Since 
the model is a free fleld theory, the couplings simply scale according to their canonical 
dimensions, and /^-functions can be deflned by using dimensionless parameters 

where k is a flducial scale (renormalisation point), giving 

0^ = K—— = -2m^ 
dK 



or 

P'^^-^r l30 = -(^^-iy (16) 

in the {r — 9) co-ordinate system. The vector fleld {3 — — (y — 1)9-^ — is plotted 

in flgures 1, 2 and 3 for dimensions 1, 2 and 3 respectively. The two dimensional case 
is plotted in polar co-ordinates in flgure 2 since it is clearly natural to make 6 periodic 
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in this instance, with period 27r. This is reflected in two dimensional conformal field 
theory where the field 'd ~ -\/27r<^ is not a primary field, but the operator : e*^'^ : is 
primary, with correlator 

<: e*^'^^^) :: e^^'^^^'^ :>f5i (17) 

\x — 



[7] . It appears, therefore, that the operator z = -^e is the most natural, not only 
from the point of view of conformal field theory but also from the geometrical point 
of view taken here. There is no obvious reason for making the fields periodic in the 
other dimensions. 



It seems natural to enquire if the vector flows plotted in the flgures are related in any 
way to the metric in equation (13). Clearly the integral curves of (16) in D = 2 are 
radial and thus geodesies of the metric (13). This is not true for D ^ 2. If a curve 
x°'{t) is parameterised by t and has tangent vector C,"- = the condition that x°'{t) 
be a geodesic is given by the geodesic equation 

C'VftC = AC (18) 



where the right hand side, A 0, allows for the possibility that t might not be affine 
[8]. 

If 13 in (16) is used in the geodesic equation (18) with the metric (13) one flnds the 
following conditions on r(t) and 0{t) 




If D 7^ 2, the only solutions, for < D < 4, are either r = (i.e. = 0) with 
A = — + or^ = with A = — y. Thus the bold lines in flgures 1 and 3 are 
geodesies. The other curves fail to satisfy the geodesic condition (18) for D ^ 2, and 
can be interpreted as being repulsed from the (/)-axis by a "force". This conclusion 
also holds for D = 4 where the only renormalisation group flows which are geodesic 
are = with A = —3 or = with A = —2. 



The 1-D Ising Model 



Clearly it is of interest to extend the above analysis beyond free held theory to the case 
of a non-trivial interacting theory. Consider therefore one of the simplest interacting 
theories - the one dimensional Ising model. 
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The one dimensional Ising model on a periodic lattice of N sites is described in [9] . 
The partition function is 



Zn = ^ exp 
{^} 



N N 



(20) 



where K = and h= with J the spin coupling and H the external magnetic field, 
(periodic boundary conditions require crjv+i = ci). ZN{K,h) can be conveniently 
expressed in terms of the transfer matrix 



V = 



++ 



V. 



+- 



^-K ^K-h 



as Zn = TrV^. 
Diagonalising V gives the eigenvalues 

A± = |cosh/t ± Vsinh^/i + e-4-^l 



Thus 



1+ 



A. 



A+ 



TV 



Zm — 

and in the limit of large N 

(1/N) InZjv = + In (^cosh/i + Vsinh^ h + e-^-^) , 
which is essentially the negative of the free energy per unit volume. 
Using equation (6) the line element is easily obtained 



-4K 



ds^ = 



X 



(sinh2/i + e-4^)^/^ 
4e"^^ cosh /i + 8 sinh^ h (cosh h + V^inh^T+e-^) 
(cosh h + Vsinh^ h + e-^^) ^ 

+ 4 sinh h dKdh + cosh h dh? 



dK^ 



(21) 



(22) 



(23) 



(24) 



It is convenient to change to a set of co-ordinates in which the metric is diagonal. To 
this end define p — e^^ sinh h and use K and p as co-ordinates. The metric is then 



ds^ = 



1 



V'(l + p2)(e4^ + p2) 



4 e^^rfK^ 



+ 



dp'' 



(25) 
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The Ricci scalar is easily calculated and is found to be 



n 



1 
2 



1 + 



(26) 



or, in terms of the physical co-ordinates K and h 

1 + 



Ve^^ sinh/i2 + 1 



(27) 



This expression was first obtained by Janyszek and Mrugaia, [10] , and the function 
is depicted in Fig. 4. Note that the curvature diverges at the critical point T = 0, 
{K oo). 

To investigate the behaviour of the renormalisation flow in this lattice model the fi- 
fmictions are replaced by a discrete recursive map, [11] . In its simplest form this 
map is obtained by asking: can one find new couplings K' and h' such that 



ZN{K',h') = ZN{K,h) ? 
{A is a normalisation factor.) 
Equation (28) is easily satisfied by demanding 



(28) 



\-K' 



„K'-h' 



= A' 



(29) 



giving the recursive formulae: 



g2ft' ^ ^2h cosh(2K + h) 



cosh(2K - h) 



(30) 



4K' ^ cosh(4JC) + cosh(2fe) 
2cosh^(/i) 

The normalisation factor A is unimportant for the present analysis. 

Note that the combination smh{h') = e^^ sinh(/i) is a renormalisation trans- 
formation invariant. Strictly speaking this transformation is not always invertible 
and so it is incorrect to refer to it as a renormalisation group transformation. The 
renormalisation flow is depicted in Fig. 5. 

Since the renormalisation transformation is a discrete map here, the /3-functions can- 
not be defined in terms of continous a vector field. However, it still makes sense to 
ask if any of the flow curves in Fig. 5 are geodesies of the metric (24). The answer is 
that only the special curve running between the two fixed points = and K = oo 
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is a geodesic, i.e. the seperatrix h = 0, the other flow hnes are "repulsed" from 

the i^T-axis and never cross it, but crowd together asymptotically in the manner of a 
caustic in geometrical optics or fluid mechanics. The geodesic distance between the 
two flxed points is 7r/2. This behaviour is suggestive of a dynamical interpretation of 
renormalisation flow. 



General Considerations 



One is led to enquire into the nature of the "force" and to investigate the possibility 
of a dynamical description of the renormalisation group flow in terms of an equation 
for the /3-functions. It will be shown that the /3-functions obey a dynamical equation 
analogous to the viscous flow of a fluid in a curved space under the influence of a 
potential, with a constraint which dictates that the kinetic energy equals the potential 
energy. 

The argument hinges on the renormalisation group equation for the two point corre- 
lators of the theory 

(32) 

If the couplings are scaled to be dimensionless, then all ^a{x) have canonical mass 
dimension D, thus the usual scaling argument gives 

= -p''d,{^aix)My)) - daP''{Mx)My)) 

- d,/3'{^a{x)^,{y)). (33) 



Integrating over all y and using translational invariance leads to 

P'dcGab + {dan Gch + [dk^ G,a = -DGab- (34) 

This equation states that /? is a conformal Killing vector for the metric Gab and it is 
straightforward to check this property using (13) and (16). That the renormalisation 
group equation can be written in terms of a Lie derivative was noted by Lassig in 
[12] and explored in more detail in [13] . A subtlety arises in D = 4 in that the integral 
over y diverges and 6*^2^2 is inflnite, due to the short distance singularity as y ^ x. 
However, if 4 — D = e is kept positive and the limit e ^ is only taken after the 
scalar curvature (14) is calculated, it appears that this is a co-ordinate singularity 
rather than a pathology in the geometry. 

If the Lie derivative on the left hand side of equation (34) is written in terms of the 
Levi-Civita connection for the metric Gab 

P'dcGab + {dan Gob + {dbn Gac = {V aP") G^b + {VbP") Gac (35) 
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and then equation (34) is contracted with one finds 



p'^VbP'' = -C'^VbU - DP" (36) 
where the "potential" is given hy U = \f3^(3'^Gi)c- The right hand side can be in- 

— * 

terpreted as a force, and the potential, U, pushes (3 away from geodesic flow unless 
VhU oc /3b . The second term on the right hand side of (36) plays the role of an 
isotropic frictional force. 

This dynamical picture is compatible with the analysis of geodesic flow for free fleld 
theories since in this case U = K'^j— + (t ~ -*-) (^'^''f'^^^] giving 

G"'dbU^[^'^' ^n_^YQ ) ^^^^ 

and 

G^'^dbU = - ^P"" for ^ = 

C'^dbU = -(^^-l^p" for r = 0, L» < 4 (38) 

so both these cases result in geodesic flow from equation (36). Note the special case 
of two dimensions where U = (l/2)r^ is the harmonic oscillator potential. 

The derivation of equation (36) relied only on the renormalisation group equation 
(32) and the fact that the composite operators, have canonical dimension D, and 
so is completely general, even for interacting field theories. However, the question of 
singularities in ($a(a;)$6(y)) as \x — y\_^ must be addressed. If {^a{x)^b{y)) ~ 
l^_^^|fc with k < then Gab = / d^x{^a{x)^biO)) is finite (assuming the correlator 

$a(3^)^&(0) fast enough for large \x\). But if A; > D, then Gab will be singular and 
must be regularised. Singular metric components do not necessarily indicate singular 
geometry and the curvature may be finite as the regulator is removed, as the example 
of free field theory in £) = 4 shows. 

In general, however, it may prove necessary to subtract counter-terms from the two 
point correlators. This has the effect of modifying equation (34), so that it reads 

P'dcGab + {daP'')Gcb + {dbP'')Gac = -DGab - Xab (39) 

where Xab is a new tensor introduced by the subtraction procedure (see for example 
[4]). The dynamical equation (36) is then modified by the introduction of a non- 
isotropic "friction" and becomes 

p^VbP"" = -G^^dbU - Dp"" - x^bP^- (40) 
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The condition for geodesic flow is now C'^dbU oc /3« and x^bP^ « P"- Equation 
(40) holds even for renormahsation schemes in which the /? functions have exphcit k 
dependence. 

In conclusion, it has been demonstrated that the metric (5) can give geodesic renor- 
mahsation group flow in free field theories, under certain circumstances. This has 
been interpreted in terms of the dynamical equation (36), analogous to fluid flow 
with friction in a curved space, under the influence of a potential with the constraint 
[/ = I I /3 p, similar to the virial theorem for a collection of harmonic oscillators in 
statistical equilibrium. For a general theory the extra non-isotropic frictional term 
in (40) may be necessary. Geodesic renormahsation flow has also been demonstrated 
in the one dimensional Ising model for the crossover between the two flxed points at 
T — oo and T — with vanishing magnetic held. 

Clearly it would be of particular interest to examine these ideas in a non-trivial 
interacting fleld theory, where it is to be expected that the more general dynamical 
equation (40) will be relevant and work is in progress on this. 

It is a pleasure to thank Denjoe O'Connor for many useful discussions Concerningl 
the renormahsation group. 
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Fig. 1: RG flow for free fleld theory in 1-D 
The bold hnes are geodesies 



^ 


^^5r ^ 


/ / 


s \ 



n 1 1 T" 

-4-2 2 

Fig. 2: RG flow for free field theory in 2-D 
All RG trajectories are geodesies 




Fig. 4: Ricci curvature for the 1-D Ising model 
The Une of the central ridge is a geodesic 



